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Abstract 

The program T0PAZ0 was developed for computing a variety of physical observables 
which are related to the e + e~ annihilation into fermion pairs and to the large an- 
gle Bhabha scattering around the Z resonance. Among them, the Z parameters or 
pseudo-observables, the de-convoluted cross sections and those dressed with QED 
radiation, and finally the forward-backward asymmetries. The calculations are per- 
formed both for a completely inclusive experimental set-up and for a realistic one, 
i.e. a set-up with cuts on the acollinearity angle, on the energy of the outgoing 
fermions or on their invariant mass and angular acceptance. The new version, 4.0, 
includes several innovative features. First of all, the most important new capabili- 
ties since previous versions are recently computed electroweak and QCD correction 
factors that are relevant at the Z resonance in the light of the present experimental 
accuracy. Among them, the effect of the next-to-leading O (am 2 ) corrections and 
those radiative corrections to the hadronic decay of the Z which provide complete 
corrections of O (aa s ) to T(Z — > qq) with q = u,d, s,c and b. Secondly, the program 
has been upgraded to cover two-fermion final states at LEP 2 energies, where some 
of the assumptions made for earlier versions are no longer valid. In particular, to 
this aim all the electroweak radiative corrections that are negligible at the Z peak, 
but relevant far from it, have been added for s-channel processes, e.g. purely weak 
boxes, next-to-leading O (a 2 ) and leading O (a 3 ) QED corrections. 



NEW VERSION SUMMARY 



Title of program: TOPAZO 4.0 
Catalogue number: 

Program obtainable from: CPC Program Library, Queen's University of 
Belfast, N. Ireland (see application form in this issue); also upon request to 
one of the authors. 

Reference in CPC for earlier versions of program: G. Montagna, O. Nicrosini, 
G. Passarino, F. Piccinini and R. Pittau, Comput. Phys. Commun. 76 (1993) 
328 (version 1.0); G. Montagna, O. Nicrosini, G. Passarino and F. Piccinini, 
Comput. Phys. Commun. 93 (1996) 120 (version 2.0). 

Catalogue number of previous version: ACNT 

The new version supersedes the previous ones 

Licensing provisions: none 

Computers on which this or another recent version has been tested: DEC- 
ALPHA 3000, HP-UX 9000; 

Operating system under which the new version has been tested: VMS, UNIX 

Installations: INFN, Sezione di Pavia, via A. Bassi 6, 27100 Pavia, and Sezione 
di Torino, via P. Giuria 1, 10125 Turin, Italy 

Programming language used: FORTRAN 77; exception to standard: use of 
REAL* 16 variables 

Memory required to execute with typical data: « 300 kbyte as evaluator of 
observables in seven energy points 

No. of bits in a word: 32 

No. of processors used: one 

Subprograms used: NAGLIB [[jj 



Keywords: e + e~ annihilation, Bhabha scattering, LEP, Z resonance, elec- 
troweak, extrapolated and realistic experimental set-up, QCD corrections, 
QED corrections, pure weak corrections, radiative corrections, Minimal Stan- 
dard Model, de-convoluted and realistic observables. 
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Nature of physical problem 

An accurate theoretical description of e + e~ annihilation processes and of 
Bhabha scattering for centre of mass energies at the Z resonance (LEP 1) 
and above (LEP 2) is necessary in order to compare theoretical cross sections 
and asymmetries with the experimental ones as measured by the LEP collabo- 
rations (realistic observables) . In particular a realistic theoretical description, 
i.e. a description in which the effects of experimental cuts, such as maximum 
acollinearity, energy or invariant mass and angular acceptance of the outgo- 
ing fermions, are taken into account, allows the comparison of the Minimal 
Standard Model predictions with experimental raw data, i.e. data corrected 
for detector efficiency but not for acceptance. The program takes into ac- 
count all the corrections, pure weak, QED and QCD, which allow for such a 
realistic theoretical description. The program offers also the possibility of com- 
puting the Z parameters (pseudo-observables) including the state-of-the-art of 
radiative corrections, which is important for the indirect determination of the 
fundamental Standard Model parameters. 

Method of solution 

Same as in the original program. A detailed description of the theoretical 
formulation and of a sample of physical results obtained can be found in || . 

Summary of revisions 

• In 1995 the CERN Report on Precision Calculations for the Z reso- 
nance |J provided as basic documentation the theoretical basis for up- 
grading existing calculations and FORTRAN programs. 

Although the '95 analysis remains quite comprehensive, an update of the 
discussion of radiative corrections has become necessary for one very good 
reason: a sizeable amount of theoretical work has appeared following 
ref. p|. In particular, a crucial amount of work has been performed in 
providing higher-order corrections. 

In ref. [|J the two- loop O (a 2 m 2 ) are incorporated in the theoretical cal- 
culation of M w and sin 2 6^. More recently the complete calculation of 
the decay rate of the Z have been made available ||. The only case 
which is not covered is the one of final b quarks, because it involves spe- 
cific O (a 2 m 2 ) vertex corrections. For simplicity the above correction 
factors will be referred as sub-leadings. Another recent development in 
the computation of radiative corrections to the hadronic decay of the Z 
is contained in two papers which together provide complete corrections 
of O (aa s ) to T(Z — > qq) with q = u,d, s,c and b. In the first reference 
of || the decay into light quarks is treated. In the second one the re- 
maining diagrams contributing to the decay into b quarks are considered 
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and thus the mixed two-loops corrections are complete. 

TOPAZO [0, [|, which was involved in the analysis of refs. [§, |{J, has been 
constantly updated and we focus, in this note, on a short description of 
the most important new capabilities since the '95 version. 

— The crucial upgrading refers to the allowed values for 0U0, a flag 
contributing to the estimate of the theoretical error. The old option 
0U0 = 'Y'/'N' has now 3 different entries. The effect of the old 
choice ' N' (no re-summation at all of the bosonic corrections, the 
extreme option) is disappeared. Instead we have now 0U0 = 'N' 
or ' L ' or ' S ' . If one wants to compare with the old default then 
0U0 = ' N ; (new default) should be selected, 0U0 = 'L' is a new 
variant of the re-summation without sub-leadings and could be used 
to estimate the uncertainty before the introduction of sub-leadings, 
finally 0U0 = ' S ' includes the implementation of the [|J correction 
factors and represents the recommended choice. TOPAZO is based an 
a series of theoretical options OUn, n = 0,. . . ,7. All these options 
(n > 0) are still active for 0U0 = ' N ' or ' L ' but most of them are 
internally de-activated for 0U0 = ' S', thus decreasing the overall 
uncertainty. The default setting remains 

0U1 =' Y', 0U2 =' N', 0U3 =' Y', 0U4 =' N', 
0U5 =' N', 0U6 =' Y', 0U7 =' N' 

Changing the parameters is governed by a FORTRAN job-control file, 
typical example is: 

CHARACTER* 1 OVAL 
CHARACTER*3 OFFS 
0FFS= '0U0' 
0VAL= 'S' 

CALL TCFLAG(0FFS , OVAL) 

— There is a major change with respect to the internal philosophy of 
TOPAZO. A new subroutine has been introduced, 

TBASIC (NT , ST , FVECMS , IFLAG) 

Once TOPAZO is called then the first step is governed by an internal 
call to subroutine TWIDHTO. Here, as a primitive calculation, both 
M w - the W boson mass - and s 2 - the improved one-loop solution 
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for the sinus of the weak mixing angle - are obtained as iterative 
solutions of the two relevant renormalization equations. After that, 
all Z-parameters are computed. As a matter of fact, most of the 
users are solely interested in deriving pseudo-observables. For this 
reason, it is now possible to call TOPAZO with OEXT = ' P ' where the 
calculation of the realistic observables is skipped with a considerable 
gain in CPU time. 

In computing M w , s 2 and the rest of the pseudo-observables there 
are essentially five functions which are touched by the inclusion of 
the sub-leadings. Each of these functions depends on M H ,m t ,s 2 
and M w . As a function of the ratio M H /m t they are given in [|| in 
terms of two expansions, M H /m t ^> 1 and M H /m t <C 1. In between 
it is up to the reader to decide what to do. In between, however, is 
exactly where the the present experimental data would like to see 
the Higgs boson. In order to avoid any kink in the evaluation of 
pseudo-observables one has to interpolate numerically. The origin 
of the kink is related to the two-loop correction factor, Ap^ 2 \ which 
in the hierarchy of the effects is the dominant one. We decided to 
interpolate with some accuracy for this function. 

The most important upgrading in the electroweak/QCD interplay 
is represented by non-factorizable QCD and EW corrections to the 
hadronic Z boson decay rate. The Born result receives both QCD 
and EW corrections and, so far, one used factorization 

r = rw fi + - 

However the correct implementation of the new result || is giving 
us 

AT(Z -> u,d,s,c) « -0.59(3) MeV, (2) 

or 

, , Ar (Z -> hadrons) 0.50 

^"■» a - r(Z^hadron 5 ) ~*m3~°-° 0L < 3 > 

There is no space enough here to account for a complete descrip- 
tion of the QCD corrections implemented in TOPAZO. We simply 
note that subroutine TCORRQCD has been completely re-designed and 
function RRUNM allows the running quark masses to be computed at 
the corresponding pole mass. For instance subroutine TCORRQCD 
computes 
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- the running of the s-quark mass up to the c- and 6-quark thresh- 
old 

- the c-quark mass at the c-quark threshold 

- the c-quark mass at the 6-quark threshold 

- and finally the running c-quark mass at any scale. 

- The 6-quark mass at 6-quark threshold 

- and finally the running 6-quark mass at any scale. 

In the following we present a list of the relevant corrections: 

1. The O (a s ) corrections including non- vanishing quark masses. 

2. Flavor non-singlet corrections in the massless limit are the same 
as those for the vector current. Whereas flavor singlet contribu- 
tions for vector currents arise only at order they are present 
already in second order for the axial part. 

3. Massless non-singlet corrections. 

4. Quadratic massive corrections. 

5. Quartic massive corrections. 

6. Power suppressed top-quark mass correction. 

7. Singlet axial corrections. 

8. Singlet vector correction. 

— It is important to discuss the most relevant variations in the predic- 
tion for pseudo-observables. We refer to the situation presented in 
refs. || For this reason we have taken again M z = 91.1888 GeV, 





NEW versus OLD 


TOPAZO 2.0 


TOPAZO 4.0 


Rel. shift (per- mil) Abs. shift 


M w [GeV] 


80.310 


80.308 


-0.03 


-2.1 [MeV] 


sin^ e l eS 


0.23200 


0.23209 


0.39 


8.9xl0~ 5 


T z [MeV] 


2497.4 


2496.1 


-0.052 


-1.29 [MeV] 



Table 1: Comparison of two TOPAZO's versions. Here M z = 91.1888 GeV, 
m t = 175 GeV, a s (M 2 z ) = 0.125 and M H = 300 GeV. 

m t = 175 GeV, a s (M 2 z ) = 0.125 and M H = 300 GeV (the '95 input 
parameter set) and compared some of the predictions in Table 1. 

— Purely weak boxes have been introduced for s-channel processes. 
Their effect has been already discussed at length in ||10| . 
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• QED corrections — For both the cross sections and forward-backward 
asymmetry, new QED corrections have been added. They include both 
corrections that are relevant at the Z peak, in view of the present ac- 
curacy of experimental data, and correction that are negligible at the Z 
peak, but become relevant for centre of mass energies above it. 



Cross sections - New corrections, already known in the literature, 
but not implemented in the previous versions, and completely new 
ones have been inserted. They are: 

* initial-state next-to-leading O (a 2 ) corrections to the radiator, 
according to eqs. (3.19) and (3.20) in |TT[; this radiator is se- 
lected by OHC = Y and ORAD = A; 

* additional initial-state higher-order corrections to the radiator, 
according to eqs. (3.29) and (3.30) in fill ; this radiator is se- 
lected by OHC = Y and ORAD = D; 

* additional initial-state higher-order corrections to the radiator, 
including part of the leading O (a 3 ) corrections, according to 
eqs. (3.31) and (3.32) in Oj; this radiator is selected by OHC = 



Y and ORAD = E; 
* the full radiator including initial-state next-to-leading O (a 2 ) 
and complete leading O (a 3 ) corrections, according to the ana- 
lytical results obtained in [fL2| ; this radiator is selected by OHC 
= Y and ORAD = F (recommended choice). 

For each of these choices, also the primitives have been computed 
analytically and implemented in FUNCTION TPRAD. Moreover, also 
the O [a 3 ) contributions to the electron structure function and its 



primitive, according to [13|, have been implemented in SUBROUTINE 



TSTRUCFUN and SUBROUTINE TPDFUN. 

Forward-backward asymmetries - For OHC = Y, additional contribu- 
tions to the forward-backward radiator, that become relevant above 
the Z peak, are included, according to eqs. (43), at the O (a), and 
(46), at the O (a 2 ), in [14]. These contributions are active only for 



s-channel asymmetries. Also in this case, the primitive has been 
computed analytically and implemented in FUNCTION TPRAD. 



OHC = N sets automatically the radiator adopted in the previous versions 
of the program, both for cross sections and asymmetries. 

The effects of the above QED corrections have been investigated in ||12|| . 
They can be summarized as follows: at LEP 1 the initial-state next-to-leading 
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O (a 2 ) corrections are negligible, whereas the leading O (a 3 ) ones introduce a 
negative shift of the total cross section of about 0.07%. At LEP 2, the two 
corrections tend to compensate one another, leaving a net effect of around 
0.2-0.4% when the Z radiative return is included and about 0.1% otherwise. 

It is worth emphasizing, that in order to implement the new radiators and 
their primitives, a new library of Nielsen's polylogarithms has been created, 

SUBROUTINE TP0LYL (X , EP , SI 1 , S12 , S13 , S21 , S22) 

It returns Si ; i(x), . . . S 2 ,2{x)- 

The running is governed by a FORTRAN job-control file, where the necessary 
input parameters are fixed. The program calls SUBROUTINE TINIT (sec below), 
where all the flags are initialized to their default value. It is allowed to change 
the default settings by calling the auxiliary SUBROUTINES TCUTSET , TCFLAG 
and TC0PT. In particular, the theoretical options OUn are re-set by subroutine 
TC0PT, kinematical cuts by subroutine TCUTSET and all the remaining TOPAZO's 
flags by subroutine TCFLAG. A change of a non-existing flag is signalled by a 
warning message: 

T0PFLAG: FLAG NOT RECOGNISED: FLAGNAME 

Then SUBROUTINE T0PAZ0 is called, according to the following calling state- 
ment: 

SUBROUTINE T0PAZ0 (NRTS , RTS , ZMT , TQMT , HMT , ALST , 0TPP0 , 0TPR0) . 

The first six entries are input parameters, whose meaning is 
NRTS : number of energies; 

RTS : array dimensioned as RTS (NRTS), containing the values of the energies; 

ZMT : the Z-boson mass (GeV); 

TQMT : the top-quark mass (GeV); 

HMT : the Higgs-boson mass (GeV); 

ALST : the value of a s {M 2 z ). 

The last two variables are output quantities: 

0TPP0 : array dimensioned as 0TPP0(24), containing the values of the pseudo- 
observables; the internal ordering is: 



hadronic peak cross-section 
partial leptonic widths 



mass of the W 



M w (l) 

Chad (16) 

r,(2),r e (3),r M (4),r r (5) 
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partial hadronic widths T u (6), 1^(7), r c (8), T s (7), Tfe(9) 

the total width T Z (U) 

the total hadronic width T/ l (19) 

the total invisible width r int) (20) 

ratios Ri(l5), R b {17), R c (22) 

asymmetries and polarization A 1 (12), A (13), A 6 (18), 



^ B (21),4 R (23),^ R (24) 
effective sinuses sin 2 6^(10), sin 2 9 b eS (11) 



OTPRO: array dimensioned as OTPRO (26*NRTS) , containing the values of the 
realistic observables and relative numerical errors: 

a e ± Aa e 0^ ± Ao^ a T ± Aa T 

Chad ± A^had 

RO { a b ± Aa 6 .Rf, ± AR b a c ± Act c i? c ± Ai? c 
A e ±AA e ±AA^ A T ±AA T 

FB FB FB FB FB FB 

A c ±AA C A b ±AA b 

FB FB FB FB 

Results are stored in array OTPRO(K) according to K = 13 (2(I-1)+J-1)+L 
with 1 = 1, NRTS running over energy points, L = 1,13 running over the type 
of realistic observable (cr e ,(x^ etc.) and J = 1,2 for the central value and the 
numerical error. 

The output is still governed by the job-control file, by means of the printing 
of the arrays 0TPP0(24) and OTPRO (26*NRTS), to be declared in COMMON. The 
output from the program is self-explanatory. 

All the internal flags of the program are initialized by means of SUBROUTINE 
TINIT, which is called by the job-control file. Here the complete list of inputs 
is given, devoting particular care in commenting the meaning of the new flags. 
For more details concerning flags already present in earlier versions, the reader 
is referred to the documentation of the previous releases. 
SE : scaling factor for numerical integration error; 

OWEAK : residual weak corrections: running (R) or fixed (F); at LEP 2 the 
correct choice is R; 

OWBOX : weak boxes included (Y) or not (N) . Away from the Z-resonance the 
correct choice is Y; 

QUO : includes (S) the sub-leading two-loop O (a 2 m^) corrections according 
to {I . It is the recommended choice; 

OHC : next-to-leading and higher order hard photon contribution included (Y) 
or not (N) ; 
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ORAD : selects the type of next-to- leading and higher order hard photon contri- 
butions to be included (A,D,E,F). F is the recommended choice; 
DAL : selects the value of a^(JI^): 128.896 (Y), 128.87 (N), user-defined (V); 
if OAL = 'V then a new value for a~^(M^) is selected through the calling 
procedure CALL TCAQED (ALPHANEW) . 
DAAS : the scale in a s in mixed aa s corrections; 

DNP : array dimensioned as ONPT(NRTS) selecting (Y) or not (N) the pair pro- 
duction correction for each energy; 

ONIF : array dimensioned as ONIF(NRTS) selecting (Y) or not (N) the effect of 
initial-finale state QED interference for each energy; 

□EXT : selects the extrapolated (E) or cut (C) branch for \i and r. With OEXT 
= 'P' pseudo-observables are computed and control is returned; 
OFB : selects (Y) or not (N) different treatment of cuts for cross sections and 
forward-backward asymmetries for /i and r; 

□CUTS : array dimensioned as OCUTST(NRTS) ; if OFB = Y selects (Y) or not (N) 
a cut on the invariant mass after initial-state radiation for \i and r cross sec- 
tions; 

DCHAN : selects the full Bhabha scattering matrix element (F) (default) or only 
the s-channel part (S); 

□CUT : array dimensioned as OCUT(NRTS); selects a cut on s': no cut (NC), 

cut on hadronic channels only (HC), cut on all channels (FC). Here s' is the 

invariant mass of the e + e~ system after initial state radiation; 

□CUTF : as OCUT, but for final-state fermions invariant mass; 

OTHRMT : if OEXT = C, selects invariant mass (M) or energy (E) threshold for /x 

and r; 

OCREE : if OEXT = E and OCHAN = S, i.e. for s-channel electrons, it allows for 
an s' cut; 

OCUTES : array dimensioned as OCUTES(NRTS); if OCREE = Y, it allows (C) or 
not (N) for an s' cut for each energy point; 

OTHRE : selects invariant mass (M) or energy (E) threshold for electrons; 

OFS : selects the treatment of higher order final-state QED corrections; default 

(D) or perturbative (Z); 

OCN: includes (Y) or not (N) the contribution of electromagnetic jets for final- 
state calorimetric electrons; 

XMED : fictitious separator for 1-dimensional integrations to improve numerical 
convergence; 

OVNAL : if OAL = V, enters the user-defined value a~^(M 2 z ) = OVNAL; 

SC: if OAAS = Y, defines the scale in a s in mixed aa s corrections; 

ZPCUT : the minimum fraction of squared invariant mass of the final state after 

radiation of the additional initial-state pair; 

OXCUTS : if OCUTS = Y, it enters the value of the s' cut, s /s; 



10 



DXCUT : if OCUT = Y, it enters the value of the s' cut, s /s; 

OXCUTF : if OCUTF = Y, it enters the cut on the value of invariant mass squared 

of the final-state fermions; 

SOCUT : array dimensioned as S0CUT(3); SOCUT(l) is the minimum invariant 
mass (GeV) of the final-state electrons, used if OTHRE = M; SOCUT (2) and 
S0CUT(3) are the same as SOCUT(l) for // and r respectively, used if OTHRMT 
= M; 

EO : array dimensioned as E0(3); E0(1) is the minimum energy (GeV) of the 
final-state electrons, used if OTHRE = E; E0(2) and E0(3) are the same as 
E0(1) for fi and r respectively, used if OTHRMT = E; 

THMIN : array dimensioned as THMIN(3); it enters the minimum scattering an- 
gle (deg) for electrons, ji and r, respectively {symmetrical angular acceptance 
is understood); THMIN (2) and THMIN (3) used only if OEXT = C; 
THMINP : the same as THMIN, for the antiparticles; 

ACOLL : array dimensioned as AC0LL(3); it enters the maximum acollinearity 
angle (deg) for electrons, \i and r, respectively; ACOLL (2) and ACOLL (3) used 
only if OEXT = C; 

OXCUTES : if OCUTES(I) = C, it enters the s /s threshold for the I-th energy 
point; 

DEL : the semi- aperture of the electromagnetic jet (deg); used only if OCN = 
'Y\ 

The internal default has been set to the following choices: 

OWEAK =' R', OWBOX =' Y', OHC =' Y', ORAD =' F' 
XMED = 0.98D0, OAL =' Y', OCN =' Y', DEL = 0.5D0 
ZPCUT = 0.8D0, OFS =' D', ONP(l) =' Y', ONIF(l) =' Y' 

A change for some of the kinematical cuts is illustrated by the following 
calling sequence: 

IND= 1 
DMY= -1.D0 
ANG= 40. DO 
ACL= 10. DO 
ETH= 5. DO 

CALL TCUTSET (IND.DMY, ANG , DMY , ACL , ETH , DMY , DMY , DMY) 

With IND = 1 only cuts for the electron will be touched and moreover DMY 
= -1, i.e. DMY < implies that variables SOCUT(l), THMINP (1), OXCUT, 
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OXCUTF and ZPCUT keep the default value, while variables THMIN(l) , ACOLL(l) 
and EO(l) are reset to new values 40°, 10° and 5 GeV. 

After initialization a call is performed to subroutine TWIDTHO after which 
the control is returned if OEXT = ' P ' , otherwise subroutine TEWEXT and TEWCUT 
are called and the evaluation of realistic observables starts. In TEWEXT weak 
corrections are computed and the physical quantities are convoluted when no 
cuts (but on the invariant mass) are applied. In TWCUT weak corrections are 
computed and the physical quantities are convoluted for leptons, when cuts 
are applied. 

The previous versions of T0PAZ0 can be used to compute Z-parameters and 
de-convoluted observables but also to obtain predictions for QED-dressed dis- 
tributions, over some realistic set-up, resembling the experimental raw data. 
The most severe limitation of those versions is that they have been developed 
by having in mind the accurate evaluation of theoretical observables at the 
Z peak. The new version of the program, T0PAZ0 4.0, offers remarkable im- 
provements in two directions. First, towards the goal of taking into account 
all those radiative corrections that can become relevant, when we consider the 
experimental accuracy reached at present by the LEP 1 experiments. Of the 
same relevance is the effort of adding those radiative corrections that are negli- 
gible for a centre-of-mass energy around the Z peak, but that become relevant 
for energies above it. 

Restrictions on the complexity of the problem 

Analytic formulas have been developed for an experimental set-up with sym- 
metrical angular acceptance. Moreover the angular acceptance of the scat- 
tered antifermion has been assumed to be larger than the one of the scattered 
fermion. The prediction for Bhabha scattering is understood to be for the 
large-angle regime. Initial-state next-to-leading O (a) QED corrections are 
treated exactly for an s' cut, in the soft photon approximation otherwise. This 
means that for center of mass energies sensibly above the Z° peak (typically in 
the LEP 1.5 - LEP 2 regime), the theoretical accuracy of the C branch is un- 
der control (theoretical error < 0.3%) when excluding the Z radiative return, 
whereas including it the theoretical error can grow up to some % depending 
on the final state selected |TJ|]. In the same energy range, large angle Bhabha 
scattering becomes a t-channel dominated process: since all the QED correc- 
tions implemented are strictly valid for s-channel processes, this means that 
large angle Bhabha scattering off the Z resonance is treated at the leading 
logarithmic level. 

Typical running time 

This depends strongly on the particular experimental set-up studied and on 
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the energy range. As evaluator of realistic observables in seven energy points 
around the Z° peak, between 10 (extrapolated set-up) and 270 (realistic set- 
up) CPU seconds for HP-UX 9000. Anyway, for the realistic observables the 
CPU time depends strongly on being at LEP 1 or LEP 2, and on the scaling 
factor SE controlling the accuracy of numerical integrations. For the evaluation 
of pseudo-observables the program runs much faster. 

Unusual features of the program 

Subroutines from the library of mathematical subprograms NAGLIB |]T| for 
the numerical integrations are used in the program. 
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work before publication of the work itself. 
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Test Run Output 

The typical calculations that can be performed with the new version of the 
program are illustrated in the following example. The FORTRAN job-control 
file computes both pseudo-observables and realistic observables with OEXT = 
'E\ For the large angle Bhabha observables we use the following selection 
criterion: 

- minimum scattering angle for the e~, | cos#~| < 0.7; 

- maximum acollinearity, 9 aco \ = 10°; 

- energy thresholds, E(e ± ) = 1 GeV. 



The realistic observables are computed for three typical energies, ^/s = 
M z , 136.22 GeV and 172.12 GeV. The input parameter set is M z = 91.1867 GeV, 
m t = 175.6 GeV, M H = 300 GeV and a s (M 2 z ) = 0.120. 



W MASS (GEV) 
NU 

ELECTRON 
MU0N 
TAU 
UP 

DOWN (STRANGE) 

CHARM 

BOTTOM 

SIN~2(E) 

SIN~2(B) 

A_FB(L) EFF. 

A_LR EFF. 

TOTAL WIDTH (GEV) 

G_H/G_E 

SIGMA0.H (NB) 

G(B)/G(HAD) 

A_FB(B) 

HADR0NIC WIDTH (GEV) = 

INVISIBLE 

A_FB(C) 

G(C)/G(HAD) 

A_LR(B) EFF. 

A_LR(C) EFF. 



0.803121E+02 
0.167123E+00 
0.839180E-01 
0.839173E-01 
0.837263E-01 
0.299863E+00 
0.382626E+00 
0.299806E+00 
0.375542E+00 
0.232068E+00 
0.233372E+00 
0.152848E-01 
0.142757E+00 
0.249338E+01 
0.207399E+02 
0.414744E+02 
0.215772E+00 
0.999354E-01 
0.174045E+01 
0.501370E+00 
0.713261E-01 
0.172258E+00 
. 934328E+00 
0.665952E+00 
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E_CM (GEV) 



0.91187E+02 



SIGMA (E) 


(NB) 




0. 


. 1016832E+01 


+/- 


0. 


.6896812E- 


-05 


SIGMA (MU) 


(NB) 


— 


0. 


. 1471494E+01 


+/- 


0. 


.4361508E- 


-07 


SIGMA (TAU) 


(NB) 




0. 


. 1468142E+01 


+/- 


0, 


.4351378E- 


-07 


SIGMA (HAD) 


(NB) 


= 


0. 


. 3034906E+02 


+/- 


0, 


.3921926E- 


-06 


OTPMA 
D.LW1A {DJ 


(Nil) 




U 


. bb4ozoorj+Ul 


+/ 


A 
U 


-1 QQQ A A -1 T7 

, 1ooo141£j 


A£ 

Ub 








A 

U , 




+/ 


A 

u . 


, do lb f o4li 


AO 

Uo 


OTP MA 

bldMA (U; 


(MB; 




A 


. bzobloob.+Ul 


I / 

+/- 


A 
(J . 


, 14ooU41h- 


AC 

-Ub 








A 


H -70/1 AVCT _i_ AA 

, 1 1 z4y f bh+uu 


X / 

+/- 


A 
U . 


. DO f lUO ( Ei" 


AO 


A_rri(.r.J 






a 
U , 


-i A7/i /i /i crxnn 
. 1U ( 444b.h+UU 




A 

u . 


, f ouyyDyti 


AC 

Ub 


A UD /'mtt^ 
A_rti U1U J 






A 

u , 


. loy j boot Uo 


x / 


A 

u . 


, yUD lOLliL 


U f 


R_f D { LRU J 






u , 


1 ^1 ohoac—ao 
. iDlzUzUr. UO 




A 
U . 


, yu 1 4oDoti 


.(Y7 
U ( 


A TTRffM 






A 

u , 


cQccoonTT— n 1 
. oyDDzzut Ul 


+/ 


A 
U . 


c -1 77n77F- 


.A7 

U 1 


A UT3 /'"D > \ 






A 
(J , 


n/i i oci at? a i 
, y41obiy.h-Ul 


x / 


A 
U , 


o-l Cyl7C AC 

, zlD4/oU£i- 


AV 

-U f 


h_01 v l ^Wiv; 








U . lobzzn,+Uo 










SIGMA (E) 


(NB) 


= 


0. 


.3910009E-01 


+/- 


0. 


.3801224E- 


-05 


SIGMA (MU) 


(NB) 




0. 


.7110528E-02 


+/- 


0. 


.9921073E- 


-07 


SIGMA (TAU) 


(NB) 


_ 


0. 


.7108396E-02 


+/- 


0. 


.9910925E- 


-07 


SIGMA (HAD) 


(NB) 


= 


0. 


.6284555E-01 


+/- 


0. 


.9266594E- 


-06 


CTPMA f 13 ^ 
DlURA. {DJ 


UNts; 




a 
U . 


- lloV zIOEj Ul 


x / 

+/ 


A 

U . 




Uo 


D 13 






A 

U , 


, loUyooyri+UU 


x / 

+/ 


A 

u , 




Uo 


O T P M A f H\ 


(jnb; 




a 
U . 


. I4z (y ( ( h-Ul 


x / 

+/- 


A 

u , 


. ooobb f lh- 




D P 






A 
U , 


. zz / ZzUlri+UU 


x / 


A 

u . 


. oo4yooiEi 


AC 

Uo 


A_Jtd {e*J 






A 

U , 


QA1 QQ/17UxAO 

. oUiyo4 ( h+UU 


x / 

+/ 


A 

u . 


. 1 ( OUODOEi 


Uo 


A_FB (MU) 









. 6970068E+00 


+/- 


o , 


. 1225792E- 


A A 

-04 


A_FB(TAU) 









.6967559E+00 


+/- 


0. 


. 1224640E- 


-04 


A_FB(C) 






0. 


.6910530E+00 


+/- 


0. 


. 2064462E- 


-04 


A_FB(B) 






0, 


.4956733E+00 


+/- 


0. 


.2387142E- 


-04 


E_CM (GEV) 








0.17212E+03 










SIGMA (E) 


(NB) 




0. 


.2439700E-01 


+/- 


0. 


.3792062E- 


-05 


SIGMA (MU) 


(NB) 




0, 


.3878666E-02 


+/- 


0, 


.2972920E- 


-07 


SIGMA (TAU) 


(NB) 




0, 


.3878186E-02 


+/- 


0, 


.2971010E- 


-07 


SIGMA (HAD) 


(NB) 




0, 


.2754651E-01 


+/- 


0, 


.2768351E- 


-06 


SIGMA (B) 


(NB) 




0, 


.4584642E-02 


+/- 





. 1347729E- 


-06 


R_B 






0. 


. 1664328E+00 


+/- 


0. 


5170563E- 


-05 
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SIGMA (C) 
R_C 

A_FB(E) 

A_FB(MU) 

A_FB(TAU) 

A_FB(C) 

A_FB(B) 



(NB) 



0.6889562E-02 +/- 
0.2501065E+00 +/- 
0.8091031E+00 +/- 
0.6097632E+00 +/- 
. 6095878E+00 +/- 
0.6717413E+00 +/- 
0.5592299E+00 +/- 



0.1053195E-06 
0.4575542E-05 
0.2815138E-03 
0.5588466E-05 
0.5584414E-05 
0.1127397E-04 
0.1854536E-04 
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